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Abstract 

We report measurements of H~ and H + hyperon absolute yields as a function 
of rapidity in 158 GeV/c Pb+Pb collisions. At midrapidity, ^ = 2.29 ±0.12 
for E~, and 0.52 ± 0.05 for E + , leading to the ratio of H + /~T=0.23±0.03. 
Inverse slope parameters fitted to the measured transverse mass spectra are 
of the order of T « 300 MeV near mid-rapidity. The estimated total yield of 
H~ particles in Pb+Pb central interactions amounts to 7.4±1.0 per collision. 
Comparison to E~ p roduction in properly scaled p+p reactions at the same 
energy reveals a dramatic enhancement (about one order of magnitude) of H~ 
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production in Pb+Pb central collisions over elementary hadron interactions. 
25.75,12.38.Mh 
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Ultra-relativistic nuclear collisions provide a unique opportunity to create and study a 
hypothetical new phase of nuclear matter: the quark-gluon plasma (QGP). For an overview 
see |TJ and references therein. In the color-deconfmed QGP phase, a g eneral enhancement 
of strange-sector particle yields with respect to that from hadronic interactions, has been 
predicted @,[|. This stems from the argument that the high parton density and lower 
energy threshold for ss quark-pair produ ction in the plasma, compared to a hadron gas, 
may lead to far shorter equilibration times in the plasma than in a hadron gas. Thus, a 
high characteristic strangeness yield would result from a plasma phase essentially unaltered 
during the subsequent hadro nic expansion. Moreover, in an environment with high baryon 
density, as e.g. in Pb+Pb collisions at the SPS 0, ss pair formation should be favored if 
the lowest available u- and <i-quark energy levels are above 2m s . 

The strangeness enhancement and its relation to QGP formation have been studied 
within the QCD framework [|5]||, confirming earlier predictions from kinetic equilibrium 
models In lowest-order perturbative QCD, ss quark pairs are created via gluonic 

fusion (gg — > ss) and light quark-antiquark annihilation (qq — > ss). In a kinetically equili- 
brated partonic gas at sufficient temperature (T pa m s ), the gluonic production dominates 
(~90%), leading to chemical equilibration times comparable to the expected plasma lifetime 
||. Relaxation times for chemical equilibrium via qq annihilation are too slow The 
gluonic degrees of freedom, transiently present during collisions, would then result in a high 
abundance of strangeness. Studies of strangeness production in central A+A collisions may 
thus help to ascertain the existence of a transient partonic phase. 

Multi-strange hyperon production is expected to be particularly sensitive to the rapid 
flavor equilibration in the early stages of the collisions, when the energy density is highest 
and the phase transition is most likely ||. Multi-strange baryons and antibaryons should 
coalesce easily during hadronization of the QGP and should survive subsequent hadronic 
interactions since their inelastic cross sections are generally small. Thus, a QGP transition 
should significantly enhance multi-strange hyperon yields. In fact, a recent publication of A, 
E, and Q particle yields near mid-rapidity by the WA97 Collaboration || has demonstrated 
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such an enhancement over p+Pb collisions. In a complementary study we obtained estimate 
s of the total An yields of and E + . We shall derive an enhancement of about a factor 
10 of the S~ particle yields in central Pb+Pb collisions over and above properly scaled pp 
collision yields at similar energy. 

Studies of S~ and S + (S=±2) production in Pb+Pb collisions at 158 GeV/c at the CERN 
SPS were carried out by the NA49 collaboration using two independent analysis techniques. 
The NA49 experiment uses four large Time Projection Chambers (TPC's) for tracking and 
momentum analysis. Two of them, Vertex TPC-1 and -2 (VTPC1,2) are placed inside 
a magnetic field (with a total bending power of 9 T-m), and two Main TPCs (MTPC's) 
are placed further downstream outside the magnetic field on either side of the beam. A 



detailed description of the apparatus can be found in |H|. In addition to the standard 
configuration (target in front of VTPC1) a set of central Pb+Pb events was taken with 
the target placed about 20 cm upstream of the MTPC entrance windows. H particles were 
studied in the VTPC2 detector (traditional approach, using a magnetic field and direct 
momentum measurements, target in front of the VTPC1) and in the MT PC detector 
(special configuration, target in front of the MTPC, no magnetic field). The lack of charge 
information in MTPC runs (special configuration) caused by the absence of a magnetic field 
did not allow a separation between S~ and S + to be made. Thus the sum of H~+H + was 
measured with the MT PCs. The cascade decay 

2~ — > A + 7r~ (where A — > p + 7r~) 

and the corresponding S + decay have a characteristic topology: the first decay vertex has 
the appearence of a broken line ("kink") followed by the V° (A) decay. A V° points back to a 
"kink" rather than to the primary even t vertex. A particles were identified by reconstructing 
their decays into a final state containing only charged particles. 

Measurements of H particles in heavy ion collisions in the presence of a magnetic field 
have been reported previously (WA85/94/97 JTTJ], NA35 |TI§, and NA36 |L_3| experiments) 



and the methodology is well established. For completeness, we briefly outline the two main 



steps of this analysis: 

• reconstruction and identification of As and As is done by measuring their charged 
decay products (V° decay). Momenta and charges of the decay products are obtained 
from the curvature of the trajectories, whereas momentum and mass of the parent 
are derived from energy /momentum conservation assuming the decay scheme, e.g. 

A — > p + 7T~. 

• the second cascade decay product (so-called "7r-batchelor" arising from the "kink" 
vertex) is found for each A (A) candidate by combining all positive (negative) tracks 
which did not point back to the target, with the rec onstructed A (A) in order to 
determine whether they could originate from a common vertex. For successful pairs 
the S~ (S + ) momentum is derived using momentum/energy conservation. The S~ 
and S + candidates are required to point back to the primary vertex. 

Our analysis of Vertex-TPC data follows this approach. 

The idea of detecting Ss in tracking detectors without the presence of a magnetic field 
is new, and therefore requires more explanation, with all verticies visible vs. TPC with 
no direct information on verticies) and significantly lighter environment (track multiplicity 
lower by one order of magnitude) - the methods developed by the UA5 collaboration were 
not applicable to the analysis of heavy ion interactions. A fairly detailed description of this 
method can be found in [14|]; here we only outline the main steps: 



reconstruction of V°-type vertices done the conventional way (by combining track 



pairs), but with straight, not curved, tracks |L5 . 



for each A candidate a search is performed to find the matching "7r-batchelor" . The 
match criterion was coplanarity of the 7r-batchelor, the A decay vertex and the inter- 
action vertex. 

roughly coplanar candidates together with the preliminary estimates of both vertex 
positions are fed into a geometrical fit routine |16| . The coplanarity constraint used in 
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the fitting routines turned out to be very efficient in eliminating the majority of the 
combinatorial background. 

• momenta of parent and daughter particles are reconstructed in each vertex from 
the decay angles (while their masses are assumed). This is done by requiring en- 
ergy/momentum conservation at each vertex, separately. 

• using the A momentum, calculated at the V° vertex, in the energy /momentum con- 
servation equations for the "kink" vertex, allowed the S mass to be calculated, rather 
than assumed. 

The advantages of this approach include the easier/higher-accuracy (straight-line) track- 
ing in the large MTPC volumes, the absence of so called "E(g>B distortions" resulting from 
an inhomogeneus magnetic field, and the proximity of the target to t he active volume. The 
final results are limited, however, to the combined (H~ + S + ) yields. 

The data presented here were obtained during two different runs, both with a central trig- 
ger. The trigger selected the most central 5 % of the total inelastic cross section (~7 barn) 
for the standard configuration (analysis in the VTPC2), and 7 % - for the special config- 
uration (analysis in MTPC). It corresponds to maximum impact parameters of 6=3.5 and 



4.0 fm, respectively, or an average number of 370 participants ||17|| . The phase-space coverage 
of the two detectors is chosen symmetrically on either side of mid-rapidity (y m ^=2.9): the 
MTPC covers y 6(1.7-2.7) andp^ >0.9 GeV/c (backward hemisphere), whereas the VTPC2 
covers forward hemisphere: y 6(3.1-4.1) and px >0.5 GeV/c. Both measurements together 
cover the bulk of the S particle yield. The data sample consists of 58K events from the 
VTPC2, and 240K events from the MTPC. Analysis of the data taken with the magnetic 
field resulted in 720 H~ and 138 S + candidates, whereas analysis of the data taken without 
magnetic field gave 2000 S~ + E + candidates. Various cuts were applied to remove the 



combinatorial background; most of them being geometrical. For details see JT8| . Besides the 
coplanarity requirement used in the MTPC analysis, important cuts in both analysis include 
the distance of closest approach between A decay products, the distance of H and A vertices 



from the interaction vertex, and the impact parameters of S decay products. Fig. [l| shows 
the invariant-mass distribution for E~ and S + , reconstructed in the VTPC2, before back- 
ground subtraction. A mass resolution of FWHM=9 MeV for H~ and FWHM=11 MeV for 
S + was obtained. The invariant mass spectrum of analyzed in the MTPC is broader 

(FWHM^30 MeV) due to the lack of direct momentum measurement . The background was 
found to be combinatoric and, depending on cuts, in the range of a few percent in VTPC2, 
and about 20% in MTPC. It was estimated by mixing As (As) reconstructed in one event 
with all negative (positive) tra cks from another event. The mixed events, processed in the 
same way as the real ones, resulted in a reconstructed "combinatoric" H (S + ) signal (and, 
in "combinatoric" S+S signal in case of MTPC), which was used, subsequently, for the 
evaluation of the magnitude of background corrections. The reconstruction efficiency was 
estimated by embedding ]19 simulated decays generated by a GEANT |16|] based Monte 



Carlo into raw data events. The events with embedded cascades were treated analogously 
to the real data (pattern recognition, track fitting, hyperon selection). The efficiency was 
found to depend on rapidity and pr\ therefore each bin of (y,pr) presented here has been 
corrected individually. The average overall tracking and reconstruction efficiency convoluted 
with th e phase space acceptance was ~0.8% for the VTPC2 and ~1% for the MTPC. Fig. |2| 
shows transverse-mass spectra for the S~ and S + analyzed in the VTPC2, and H~ + S + 
analyzed in the MTPC. The shape of the distributions is approximately exponential over 
the full range. For comparison (open circles) we also present the H~ + S + sum from the 
VTPC2 data sample. Both H~ + S + spectra are in agreement. The distribution is fitted 
with an exponential function in m^: 

= C(T)m T e~^ 

dniTdy 

with T as fit parameter. The normalization constant C(T) is constrained to the experimental 
yield. The inverse-slope parameters for 5 particles in each data set are similar and in the 
vicinity of 300 MeV (see Fig. |3]a,b) . This agrees well w ith the dependence of the inverse 
slope T on the particle mass in Pb+Pb collisions established by NA49 measurements [20 
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A similar dependence is reported by the NA44 collaboration The slope parameters 

increase monotonically with partic le mass (in the case of NA49: from pions to deuterons), 
suggesting the presence of collective transverse flow in nuclear collisions at SPS energies |22fl . 
Our measurements of H inverse slope parameters are consistent with the slopes reported by 
t he WA97 experiment |J for a more relaxed trigger. 

The rapidity density distribution was obtained by integrating the transverse-mass spectra 
in two rapidity bins. They were chosen symmetrically with respect to y m id=2.9. In the 
MTPC - the first bin covered y G (1.7-2.2), the second - y G (2.2-2.7) ; in the VTPC2 - the 
first bin covered y G (3. 1-3.6), the second - y G (3.6-4.1). The inverse slope parameters in both 
rapidity bins are rather similar as illustrated in Fig. |3|a,b. The left panel shows inverse slope 
parameters for H~ + H + from the MTPC and VTPC2 and the right one for H~ from the 
VTPC2 (the analogous plot for S + lacks in sufficient statistics). Full symbols are measured 
values, open symbols have been reflected at mid-rapidity. The transverse mass distribution 
was extrapolated to regions outside the NA49 acceptance by us ing the parametrization given 
above. In Fig. ^c,d the rapidity distributions for H~ + E + from the MTPC and VTPC2 
(left) and for S~ from the VTPC2 (right) are presented. Full and open symbols represent 
measured and reflected values, as before. The two indepe ndent measurements of S particles 
in the MTPC and VTPC2 agree within errors. At mid-rapidity, the + S + rapidity density 
is measured to be about 2.7, with a major fraction of this number (~2.3) attributed to 2~ 
rapidit y density (right panel). The former value appears to be ~30% higher than the WA97 
measurement [[|. The ratio at mid-rapidity is 0.23±0.03. This agrees within errors 



with the result of WA97 for somewhat less central Pb+Pb collisions ||11|| . The systematic 
errors were studied in both detectors. The stability of the dependence of final results on the 
cuts applied showed that the systematic uncertainty in the integrated yields in both cases 
did not exceed 15%. This is confirmed by the differen ces in the S~ + S + mesurements in 
the MTPC and VTPC2 presented on Fig. |3]c. The mid-rapidity multiplicities per unit of 
rapidity, dN/dy, are 2.29±0.12 and 0.52±0.05 for E~ and S + , respectively. Quoted errors 
are statistical only. 
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Making use of reflection symmetry we have systematically investigated [2^] Gaussian 
extrapolations of the rapidity distributions in Fig. [3|c,d. The width was assumed to be 



1.1. This was motivated by the 5 rapidity distribution provided by the UrQMD model [25 
and the vector meson (f> width measured in our experiment (a =l.O7) J23J]. The resulting 
estimates of total Air multiplicities are 8.2±1.1 for S" + E and 7 .4±1.0 for E~ . These 
values are not corrected for the hyperon Q~ feed-down: Q~ — > E~ + it (B.R. = 8.3%); 
however, based on the Q yields in Pb+Pb collisions reported by WA97 || , we estimate that 
this contribution is very small (less than 0.5%). The systematic error of the extrapolation 
procedure, not included in the errors quoted above, was estimated to be less that 14% by 
comparing the results obtained using different assumptions about the shape of the rapidity 



distribution [23]. As our central trigger (5% of the total inelastic cross section) corresponds 
to a mean number of 370 participating nucleons we thus obtain an estimate of 0.020±0.002 
E~s per participating nucleon and thus 0.040±0.004 per NN participant pair. In order now 
to check the degree of strangeness enhancement implied by these results we require similar 
data for elementary pp collisions. To our knowledge no such data exist in the vicinity of our 



y/s=18 GeV but at ^=6 |2§ and 63 G eV |2"7f. The E~/A ratio is 0.016±0.011 and about 



0.06±0.02Q respectively for these two energies. Assuming that this ratio increases faster 
with y/s at the lower of the two energies we do not interpolate linearly but estimate the 
2~/A ratio to be 0.03±0.01 at y/s=18 GeV. Taking the A multiplicity at this energy t o be 



0.10±0.03 from ref. we finally arrive at an estimate of a E~ multiplicity of 0.003±0.0015 
in pp collisions at \/s=18 GeV. The E~ yield per nucleon pair in central Pb+Pb collisions 
is thus enhanced by about one orde r of magnitude, over and above the yield in elementary 
collisions. A more rigorous determination of the enhancement factor will be derived from 
analysis of pp collisions at y/s=18 GeV currently in progress by NA49. 



1 At the ISR energy (-^=63 GeV) we used the measured S + /A ratio |27j as an estimate of the 
3~ /A ratio at mid-rapidity. This provides an upper limit of the H~/A ratio in 4-7T acceptance. 
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In summary, we have presented data on the absolute yields of S~ and H as a function 
of rapidity near mid-rapidity in central Pb+Pb collisions at the SPS using two independent 
analysis techniques, one of which is a novel method incorp orating topological identification 
without a magnetic field. We have reported a strong increase in the production of S~ and 
S + with respect to pp interactions. The inverse slope parameters of the transverse mass 
distribution of H particles are consistent with previously reported systematics [20f reflecting 
the presence of transverse collective hadronic expansion flow. 
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FIG. 1. Invariant mass spectra for E and S 
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FIG. 2. Transverse mass spectra for H~ and H + analyzed in VTPC2, and + H + analyzed 
in MTPC (offset by factor 10) 
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FIG. 3. a,b:Inverse slope parameters of 3~ + 3 and 3~ (right) as a function of 5 rapidity. 
c,d:Rapidity distributions of 3~ + 5 + and 3~ (right). The open symbols are the measured points 
reflected at y m id=2-9. See text for details. 
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